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We analyze pulsar fluxes at 1400 MHz (S1400) and distances (d) extracted from the Parkes Multi- 
beam Survey. Under the assumption that distribution of pulsar luminosities is distance-independent, 
we find that either (a) pulsar fluxes diminish with distance according to a non-standard power law, 
due, we suggest, to the presence of a component with Si 400 oc 1/d, or (b) that there are very sig- 
nificant (i.e. order of magnitude) errors in the dispersion-measure method for estimating pulsar 
distances. The former conclusion (a) supports a model for pulsar emission that has also successfully 
explained the frequency spectrum of the Crab and 8 other pulsars over 16 orders of magnitude of 
frequency, whilst alternative (b) would necessitate a radical re-evaluation of both the dispersion- 
measure method and current ideas about the distribution of free electrons within our Galaxy. 



PACS numbers: 97.60.Gb,95.30.Gv,98.70.Dk,41.60.-m 

Recently, superluminal polarization currents, whose 
distribution patterns move faster than light in vacuo, 
have been invoked as sources of pulsar emission [T]. This 
idea is derived from the work of Bolotovskii, Ginzburg 
and others, who showed both that such superluminal po- 
larization currents do not violate Special Relativity (since 
the oppositely-charged particles that make them move 
subluminally) and that they form a bona-fide source term 
in Maxwell's equations [U HI [S] . The validity of these 
ideas has been demonstrated in a variety of laboratory 
experiments [H |7j [9j [10] . Moreover, by extending the 
approach to a superluminal polarization current whose 
distribution pattern follows a circular orbit, it was pos- 
sible to explain quantitatively several observables from 
the Crab pulsar, including the spacing and widths of 
the emission bands at frequencies around 8 GHz, the 
maximum of the radiation spectrum, and the overall 
continuum spectrum across 16 orders of magnitude in 
frequency [Ij. Subsequently, successful quantitative fits 
were carried out for 8 other pulsars EQ and a related su- 
perluminal model reproduced the general form of pulsar 
Stokes parameters [12] . 

In this Letter, we demonstrate a further prediction for 
rotating superluminal sources; that there is a component 
of the emission whose flux S decays with distance d as 
S oc 1/d E3 Q3], rather than the conventional inverse 
square law (S oc 1/d 2 ). Our demonstration employs 
a Maximum Likelihood Method (MLM) [TS] analysis of 
pulsar observations [16]. The MLM is carried out to cir- 
cumvent the significant Malmquist bias [17] due to the 
increasing non-detection of weaker pulsars as d increases. 



The pulsar emission component with S oc 1/d is due 
to a general property of sources that travel faster than 
their emitted waves; the relationship between reception 
time and retarded time is not monotonic and one-to- 
one [5] [18]. Multiple retarded times [19], or, if the source 
accelerates, extended periods of retarded time [IS], can 
contribute to the waves received instantaneously, result- 
ing in temporal focusing, i.e., concentration of the energy 
carried by the waves in the time domain J5] EH E3 ■ This 
effect is well known in acoustics [HI E3 ED] • It is the tem- 
poral focusing from the parts of the source that approach 
the observer at the wave speed and with zero acceleration 
that leads to the S oc 1/d flux component E3E4]- Note 
that this mechanism does not violate conservation of en- 
ergy since the enhanced flux detected in some places is 
compensated exactly by diluted fluxes elsewhere [13] El ■ 
Moreover, we emphasize that the emission discussed in 
this paper arises from true superluminal motion; electro- 
magnetic disturbances (polarization currents) that travel 
faster than the speed of light in vacuo, c [3] ■ This should 
not be confused with the apparent superluminal motion 
of certain radio sources that is thought to arise from rel- 
ativistic aberration [21] . 

Our analysis of S1400 versus d assumes that the lu- 
minosity function of pulsars is uniform throughout the 
Galaxy, i.e. that the populations of pulsars at various 
distances from the Earth are similar, each containing 
a representative spread of pulsar types, sizes and en- 
ergies. We extract 1400 MHz fluxes (iSuoo) and dis- 
persion measures from the ATNF Pulsar Catalogue [16] 
(http://www.atnf.csiro.au/research/pulsar/psrcat). To 
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FIG. 1: (a) The 1109 Galactic pulsars in the Parkes Multi- 
beam Survey plotted as S1400 versus distance d, where d is de- 
termined from the NE2001 interpretation of dispersion mea- 
sure [23]- Pulsars with periods Po < 0.1 s are shown as hol- 
low points, and those with Po > 0.1 s are displayed as filled 
points. The apparent differences between the distributions 
of the two sets of pulsars may reflect the fall-off in sensitiv- 
ity of the Parkes instrument for faster pulsars (see Fig. 2 of 
Ref. |26|V (b) Cumulative population distribution in 5i4oo 
for 9 distance bins. The mean distance of each bin is given in 
the inset key. 

eliminate statistical biases from different instruments, we 
restrict the sample to the 1109 galactic pulsars detected 
using a single instrument, the Parkes Multi-beam Sur- 
vey [22] • We use the so-called NE2001 [23] model to eval- 
uate d values from the dispersion measures given in the 
ATNF catalogue; this was shown [23J [21] to give pulsar 
positions that are more consistent with the known distri- 
butions of matter in the Galaxy than previous models. 

We first show that the Parkes observations show a 
strong Malmquist bias due to instrumental issues; con- 
sequently, the MLM [15] is essential in making quanti- 
tative conclusions about the flux-distance relationship. 
Fig. 1(a) plots the Galactic pulsars from the Parkes 
Multibeam Survey as log 10 (>Si4oo) versus log 10 (d). It is 
already obvious that data are very sparse for 61400 < 
0.1 mjy. To assess whether this is an instrumental arte- 
fact, or a fundamental property of the pulsar population, 
we group the pulsar data in bins covering certain distance 
ranges (e.g. 6.0 < d < 7.0 kpc) and plot the cumulative 
distribution functions N(Si4oo) of each bin in Fig. 1(b). 
The d bins are chosen so that they cover a reasonably 
small range of d but contain a large enough population 
for meaningful statistics (~ 100 pulsars). 

Note first that all the cumulative distribution functions 
in Fig. 1(b) tend to zero at roughly the same £1400 ■ This 
strongly suggests that low-flux part of each cumulative 
distribution is representative of the roll-off in sensitivity 
of the instrument, rather than an intrinsic property of 



each pulsar population. On the other hand, the high- 
flux sides of the curves in Fig. 1(b) are likely to be more 
representative of intrinsic properties of the pulsar pop- 
ulations. As such, they should move to lower fluxes as 
d increases. This does indeed happen, but at a slower 
rate than the inverse-square law; the 75% points of the 
functions spread over roughly a factor 2.8 in S1400, even 
though the distance varies by a factor of around 3.3. This 
is a much smaller spread than that expected for the in- 
verse square law (~ 10). 

Simple analysis thus far has suggested that the Parkes 
Multibeam Survey is subject to a substantial Malmquist 
bias because it misses a large fraction of pulsars with 
$14,00 < 0.4 mjy and cuts off completely for S1400 <^ 
0.1 mjy. Both of the latter figures are of the same or- 
der of magnitude as the calculated minimum detectable 
flux of the Parkes instrument (~ 0.15 mjy; see Fig. 2 
of Ref. [53]). To make further progress, we require a 
method that attempts to compensate for missing data, 
caused by instrumental sensitivity problems, in a sys- 
tematic way. Originally, Efstathiou et al. [15] . and sub- 
sequently a number of authors, have demonstrated that 
the MLM is very suitable for such problems by applying 
it to red-shifts of very distant objects, a data set which 
is incomplete due to instrumental problems somewhat 
analogous to those of the Parkes survey [15] . The Parkes 
database is especially suitable for treatment, since S1400 
and d values are essentially independently derived. 

Our implementation of the MLM determines the prob- 
able luminosity function, 4>(L), based on the (incomplete) 
observed data, where L is the luminosity. The technique 
fits a quasi-continuous <fi(L) to the observations under 
the assumption of an instrumental cut-off |15l . Here, the 
MLM is used to determine the most likely value of the 
exponent n in the relationship S1400 oc d~ n . The in- 
trinsic luminosity function is therefore calculated using 
the 51400 and d values from Parkes Survey for each of 
the trial exponents (n tr iai € {0.5,1.0,1.5,2.0,2.5,3.0}); 
this feeds into the self-consistent determination of 4>(L) 
and the cutoff through the relationship Li = Sidf, where 
Si = S*i4oo f° r the i th pulsar at distance di. The proba- 
bility that pulsar i is observed in a flux-limited survey is 
given by 



Pi oc (j)(Li)/ 



<f>(L)dL. 



where L m i n (di) is the minimum luminosity that a pulsar 
at distance di can have to be detected. We define a like- 
lihood function C = YiiPi- Following Ref. [15], we use 
an approach that does not assume a simple functional for 
</>(£). Instead, the luminosity function is parameterized 
as N b steps: (f>(L) = cj> k , for L k - ^ < L < L k + 
with k = l,...,Nb- The maximum likelihood function 



3 



assumes the form 



A, 



ln£ = W ( L i ~ £fc) ln < 



jALH[Lj - L min (di)] > + const. (1) 



where iV p is the total number of pulsars in the Parkes 
survey. Here, W(x) = for -AL/2 <x< AL/2, and 1 
otherwise, and 

f 0, x < -AL/2 

= < (x/AL+1/2), -AL/2 < x < AL/2 
[ 1, x>AL/2. 

The parameters <pk determining the luminosity function 
are given by the self-consistent set of equations 



,AL 



H[L k 



with k = 1 , , Nf,. The above equations are solved 

iteratively to obtain the luminosity function, with the 
goodness of fit being parameterized by the relative con- 
vergence error e = — 0i-i(^)) 2 - This is basi- 
cally the relative mismatch between successive iterations; 
the smaller the value of e, the better the representation 
of the data. In all cases, the 10 — 15% or so of pulsars 
with very high intrinsic luminosity were excluded from 
the analysis to ensure a quasi-continuous distribution of 
the luminosity function. 

When the complete Parkes data set is used (Fig. 2(a), 
solid points), we find that the derived luminosity func- 
tion converges very rapidly for n tr i a i = 1.0 and 1.5 with a 
small e. The convergence to a putative luminosity func- 
tion is considerably (~ 10 5 ) worse when one assumes an 
unphysical ritrial = 0.5, 2.5 and 3.0, as well as the com- 
monly accepted inverse-square law (ritrial = 2). The rel- 
ative convergence for these values of ritrial can be some- 
what improved by restricting the analysis to a smaller 
set of pulsars, but is still not comparable to those ob- 
tained for ritrial = 1-0 and 1.5. Overall, the best com- 
bination is ritrial = 1 with 983 pulsars fitted, implying 
>Si4oo 1/d; the error with ritrial = 1.5 is somewhat 
larger, with 980 pulsars fitted. However, the good con- 
vergence for both exponents suggests that the observed 
flux may be a mixture of the S1400 oz 1/d component and 
spherically-decaying radiation, both of which are to be 
expected from a supcrluminally-rotating source [13j . In 
slower pulsars, whose light cylinders lie further away from 
the central neutron star, the superluminally-rotating part 
of the current distribution may not be dense enough to 
give rise to a dominant nonspherically-decaying compo- 
nent of the radiation pQ. 



The main assumption of our analysis thus far is that 
pulsar populations are similar throughout the Galaxy. 
This is potentially open to question if there are two dis- 
tinct populations of pulsars, especially if some property 
of each population results in different instrumental cut- 
offs. As shown in Fig. 2 of Ref. [53], the sensitivity of 
the Parkes Instrument is limited for pulsars with periods 
Po < 0.1 s; i.e., distant millisecond pulsars are harder 
to detect than equivalent slower pulsars (see Fig. 1(a)). 
This might result in an apparent d-dependent evolution 
of the characteristics of the detected pulsar population. 
Second, though recent work [TJ [TT] suggests that all pul- 
sars possess the same emission mechanism, some opine 
that millisecond pulsars form a distinct population |27j 
and might therefore possess a different luminosity func- 
tion. Both of these concerns can be addressed by ex- 
cluding pulsars with periods Po < 0.1 ms from the MLM 
fit. The hollow points in Fig. 2(a) show the result; the 
minimum e is obtained with n tr i a i = 1.5. For compari- 
son, Fig. 2(b) shows the result of running the MLM on 
only the 43 pulsars with Pq < 0.1 s; once again, the fit is 
best for ritrial — 1-5. Therefore, in spite of the separation 
of pulsars with P < 0.1 s and P > 0.1 s, the MLM 
never obtains a minimum e for ritrial = 2; the errors for 
ntriai = 1,1.5 are smaller, often dramatically so. This 
suggests that the violation of the inverse-square law by 
the pulsar population is a robust phenomenon. 

Finally, the MLM was tested on two synthetic galax- 
ies of pulsars in which 5*1400 was constrained to decay as 
d^ 1 - 5 and d~ 2 (Fig. 2(c)). As all pulsar population mod- 
els in the literature (e.g. Refs. [U1[5S]) are both contam- 
inated with assumptions involving the inverse-square law 
and involve ~ 10 — 20 adjustable parameters, we derived 
synthetic pulsar distributions that are consistent with the 
Parkes database using Bayesian methods |28 |. As with 
the Parkes data, the luminosity function was determined 
using a range of trial exponents (n tr i a i = 0.5, 1, 1.5, 2 and 
2.5). In each case, the MLM located the correct value of 
n, with an accuracy better than ±0.5 (Fig. 2(c)). This 
gives confidence in our assertion that pulsar flux data in 
the Parkes survey violate the inverse square law. 

In the various implementations of the MLM, the in- 
ferred pulsar luminosity function (e.g., Fig. 2(d)) always 
decreases monotonically [59]. As discussed above, the 
effects of instrumental insensitivity will lead to an appar- 
ent luminosity function that falls off at low Suoo- This 
problem may have led to the commonly-held belief that 
the pulsar luminosity function is a log-normal distribu- 
tion [27]. Once one compensates for the loss of data 
(Fig. 2(d)), the maximum characteristic of a log- normal 
distribution will be absent. 

The MLM [15] therefore finds that the observed 
1400 MHz flux of pulsars does not fall off as the 
conventionally-assumed 1/d , instead returning expo- 
nents of either 1 or 1.5. There are two possible conclu- 
sions; either the dispersion measure estimates of pulsar 
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FIG. 2: (a) Blue solid points: relative convergence error 
e from the MLM applied to 983 (n tri ai = 0.5, 1) or 980 
(itriai = 1.5,2,2.5,3) pulsars from the Parkes Survey versus 
trial exponent ntriai. Red, hollow diamonds: e from the MLM 
applied to 941 (ntriai = 0.5,1) or 938 (n tr iai = 1.5,2,2.5,3) 
pulsars with periods Pq > 0.1 s. (b) MLM fit for the 43 Parkes 
pulsars with Po < 0.1 s. (c) Relative convergence errors for 
synthetic galaxies in which 5*1400 oc d~ (red, hollow points) 
and S1400 oc d~ 2 (green, filled points) versus ntriai- (d) The 
inferred luminosity function for 983 pulsars with ntriai — 1; 
large circles are raw Parkes data; small diamonds are the fit- 
ted luminosity function. 

distance |23j are radically (and consistently) incorrect by 
factors of order 10, or pulsars do, in fact, have a flux that 
falls off more slowly with distance. The former conclusion 
would call into question the widely-used [H] (but admit- 
tedly flawed) NE2001 method [53] for estimating the dis- 
tances of radio sources from dispersion measure, and also 
be at variance with the presently-accepted distribution 
of matter within the Galactic plane [22] [23] . Instead, we 
suggest that S1400 for pulsars falls off more slowly with 
distance than the inverse-square law, due to the presence 
of a component whose flux varies as 1 /d. This is in agree- 
ment with the superluminal emission model of pulsars, a 
concept that has also correctly predicted the emission 
spectrum of the Crab and eight other pulsars [TJ QT] and 
reproduced other salient features (e.g apparent bright- 
ness temperatures, Stokes parameters, small extent of 
the emitting region etc.) seen in pulsar observations [15] . 
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